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The  present  work  extends  boundary  integral  equation 
techniques/  developed  in  an  earlier  paper  for  the  Dirichlet 
and  Neumann  problems  [ 6  ] ,  to  the  boundary  value  problem  of  the 
third  kind  (also  referred  to  as  the  Robin  problem)  associated 
with  the  Helmholtz  equation  in  exterior  domains.  While  the 
existence  and  uniqueness  of  solutions  of  the  exterior  Robin 
problem  have  been  established  previously  (see  e.g.  Leis  [9  ]; 
Danilova  [ 2  1 ) /  the  application  of  the  Helmholtz  representation 
and  the  subsequent  derivation  of  boundary  integral  equations  does 
not  appear  to  have  been  studied  until  now. 

-The  considerations  below  will  show  that  the  familiar  diffi¬ 
culties  at  interior  eigenvalues  arise  in  the  case  of  the  Robin 
boundary  conditions.  The  non-uniqueness  of  solutions  of  the 
boundary  integral  equation  at  interior  eigenvalues  is,  of  course, 
an  artifact  of  the  method  chosen  to  solve  the  original  boundary 
value  problem.  As  explained  in  [ 6  ] ,  supplementary  conditions 
of  some  form  are  needed.  Either  the  fundamental  solution  may  be 
modified,  as  was  done  by  Danilova  ( 2  ]  in  her  layer  theoretic 
treatment  of  the  third  boundary  value  problem  with  an  oblique 
derivative,  or  an  additional  equation  is  needed  if  the  funda¬ 
mental  solution  is  unaltered.  In  the  present  work  we  employ  the 
Green's  theorem  approach  to  boundary  integral  equations  with 
an  unmodified  fundamental  solution  and  show  that,  unlike  the 
situation  for  the  exterior  Dirichlet  and  Neumann  problems,  the 
supplementary  condition  appears  as  an  operator  equation  of  the 
second  kind  rather  than  one  of  the  first  kind.. 
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Rather  than  repeat  the  bibliographic  details  reviewed  by 
Kleinman  and  Roach  [ 6  ] ,  we  refer  the  reader  to  the  introduction 
of  that  paper.  Moreover,  we  will  adopt  the  notational  conventions 
of  that  paper.  Section  2,  below,  briefly  reviews  those  conventions 
as  they  are  appropriate  to  the  present  considerations.  However, 
the  boundary  value  problem  discussed  here,  in  contrast  to  those 
in  [  6  ]  and,  in  the  papers  of  Leis  [5]  and  Danilova  [2],  is 
generalized  so  that  the  data  and  the  multiplier  which  appears  in 
the  boundary  condition  may  be  specified  in  LJ3D)  rather  than 
Cq  (3D)  where  3D  is  the  boundary  of  the  exterior  domain. 

This  more  general  setting  requires  some  extensions  of 
theorems  which  appear  in  the  literature,  in  particular  extensions 
of  a  regularity  theorem  appearing  in  Mikhlin  [10  ]  and  of  a 
uniqueness  theorem  due  to  Leis  [  9 ] .  More  significantly,  it 
requires  a  careful  formulation  of  a  trace  theorem  suitable  for 
functions  described  by  single  and  double  layer  potentials.  Standard 
trace  theorems,  e.g.  Miranda  [11]  were  not  sufficient  to  ensure 
existence  of  the  normal  derivatives  of  double  layer  potentials 
required  for  application  in  the  present  setting.  These  matters  are 
discussed  carefully  in  section  3  and  in  the  appendix  which  contains 
some  computational  details  collected  there  so  as  not  to  interrupt 
the  flow  of  the  main  argument. 

Having  discussed  these  background  matters,  the  main  results 
concerning  the  boundary  integral  equation  formulation  of  the  Robin 
problem  are  presented  in  section  4  and  the  equivalence  of  this 
formulation  with  the  original  boundary  value  problem  is  proven. 
Section  4  concludes  with  a  result  pertaining  to  the  scattering 
problem  for  imperfectly  reflecting  surfaces. 
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2.  Formulation  of  boundary  integral  equations . 

Let  D_  denote  a  bounded  domain  in  with  boundary  3D  which 

is  a  closed  Lyapunov  surface  with  index  1.  We  denote  the  exterior 

of  3D  (i.e.  the  complement  of  D._u3D  in  R^)  by  D+. 

Let  R=R(P,Q)  denote  the  distance  between  two  typical  points 

P  and  Q  in  R^ .  A  fundamental  solution  of  the  Helmholtz  equation 
2 

(A+k  )u**0,  lmk>0,  is  a  two-point  function  of  position,  y(P,Q), 
which  for  convenience  we  write  in  the  form 

(2.1)  Y(P,Q)  :*  -e  iJcR/27rR> 


Complex  conjugates  will  be  denoted  by  a  bar,  thus 


(2.2) 


Y(P,Q)  :*  -e 


-ikR, 


2ttR. 


we  denote  by  3/3np  differentiation  in  the  direction  of  the  unit 
vector  np  normal  to  3D  at  the  point  pe3D.  A  point  PeD+  is  assumed 
to  have  spherical  polar  coordinates  (r  ,  0  ,  6)  relative  to  a 

r  r  r 

Cartesian  coordinate  system  erected  with  origin  in  D_.  We 

/\ 

emphasize  that  throughout  wu  shall  assume  that  np  is  the  outward 

A 

drawn  normal  with  respect  to  D_,  that  is  np  is  directed  from  3D 
into  D+.  Further,  we  shall  write  3/3n~  and  3/3n+  to  denote  the  normal 
derivative  when  P-*-pe3D  from  D_  and  D+  respectively. 

Let  (Sy)(P)  and  (Dy)  CP)  denote  the  effect  at  the  point 
Pel3  of  single  layer  and  double  layer  distributions  respectively 


on  3D  of  density  u  as  follows 


r  dis 

r\ 


n  •  r  r.  ,  . 

o  V  v 

**«  u  •  , 

w  r  *  v.  J 

®  V)  O  £  *>  .  I  4.,  c 

®  M  »•«  c.'  '*i  ’  fi  > 

H  ^  i  >•  *h  ^  t) 

^  C*  £3  »->  I  CL  Q  ^4 

-  (  \  n 
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(2.3)  (Su)  CP)  :*  [  u  (q)  Y  (P,q)  dS  ,  PeR3 

)  3D  q 


and 


(2.4) 


(Du)  (P) 


3D 


U  (q) 


3  Y 

*5 


(P,q)dSq,  PER 


Further,  for  pe3D  we  define 

(2.5)  (Ku)(p)5=  [  u(q)  (p,q)  dS  . 

' 3D  3np  q 

Notice  that  K  is  a  completely  continous  operator  on  L2(3D) 

(Mikhlin  [10]).  Jump  conditions  for  the  single  and  double  layer 
potentials  on  3D  can  be  expressed  in  terms  of  K  and  its  L2 ( 3D) 
adjoint  K*  which  is  given  by 

(2.6)  (K*u) (p)  *  [  y(q)  (P/q)  dS  . 

J  3D  3nq  q 

These  jump  conditions  as  well  as  continuity  properties  of  the 
single  and  double  layers  depend  on  the  smoothness  of  the  density 
U  and  will  be  discussed  in  Section  3. 

Finally,  we  remark  that  the  single  and  double  layer  operators 
with  PeD_,  can  be  considered  as  operators  with  range  in  L2(D_) 
and,  as  such,  are  continuous  operators  from  L2(3D)  to  L2(D_) 

(see  e.g.  Miranda  [11]). 
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With  this  notation,  representations  of  solutions  of  the 
Helmholtz  equation  obtained  by  applying  Green's  Theorem  or  the 
Helmholtz  representation  have  the  following  forms.  For  radiating 
wave  functions  <f>+  satisfying 


(2.7) 


(A+k2)<j>+(P)  =  0, 


PeD+, 


(P)-ik<fr+(P) 


0, 


we  obtain 


(2.8) 


-♦+  (q) 


|ir  (p'q)}  dSq 


SIrrJ  (P)  “  (D<(>+)  (p) 


2<f>+  (P) ,  PeD+, 
<J>+  (P) »  pe3D, 
0,  PeD^ • 


whilst  for  any  solution  of  the  Helmholtz  equation  in  D_,  we 
have 


{♦-  (P'q)  “  &T  <«)Y(P,q)}  dS 

3D  q  q  J 


(2.9) 


3$ 

=  (D<J>_)  (P)  -  (S^)  (F) 


0, 


=  i 


<*>_(p)  f 


2*  (P), 


P£D+, 
pe3D, 
PeD  . 
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We  will  be  concerned  with  the  following  boundary  value 
problem: 

Exterior  Robin  Problem: 


(2.10)  (V2+k2)  u+(P)  »  0,  PeD+ 


(2.11) 

(2.12) 


3u 

+  o(P)  u+(P)  =  g+(P)  ,  pe3D 


lim  r 
r  -*■<*>  P 
P 


,3u+ 


(P) 


1 


iku+ (P) ' 


=  0 


where  g+  eL^ (3D)  and  aeL^OD).  In  the  exterior  scattering  problem 

'N-.l  f  J 

q+  is  replaced  by  -  - au  where  u  is  the  known  incident  field 

Proceeding  formally  to  obtain  boundary  integral  equations 
whose  solutions  will  subsequently  be  used  to  construct  solutions 
to  the  boundary  value  problem  we  apply  the  Robin  condition  to 
the  representation  formula  (2.8)  which  we  write  as 


f 


fs5u+i 

2u+(P)  , 

PeD+, 

(2.13) 

(p)  -  (Du+) P  =  < 

u+(p)  f 

Pe3D, 

0 

. 

PeD_. 

obtain 

(2.14) 

(I+K*+Sa)  u+(P)  =  (Sg+)  (P)  ,  pe 3D 

where  we  have  also  used  the  fact  that  on  3D 
Du+  *  f?*u+. 


? 


If  this  boundary  integral  equation  has  a  solution  “+£1*2  ( 3D)  ,  then 
equation  (2.13)  may  be  employed  to  define  u+  in  the  exterior 
domain  D+  as 


(2.15)  u+(p)  =  *  S(g+-ou+)(P)  -  *  (du+) (P) ,  PeD+. 


This  function  is  clearly  a  solution  of  the  Helmholtz  equation  in 
D+.  To  insure  that  this  solution  in  fact  satisfies  the  boundary 
conditions,  we  require  that  the  function  u+  satisfy  an  additional 
relation  which  we  derive  from  the  Helmholtz  representation  by 
computing  normal  derivatives  as  follows : 


(2.16) 


lim 

P-*-p 


3u+ 

S(^ 


D(u+)  ]  (P) 


0 


or,  using  the  relation  (3.16)  below. 


.  .  3U+  3U+  .3 

(2.17)  “"JJn-  "*■  K  ( »  ■ )  (p)  —  lim_  ■wjj —  D(u^)  (p)  —  0. 

P  q  P-*-p  p 


By  utilizing  the  boundary  condition  we  have 


(2.18)  (-g+  +  au+  +  K(g+  -  a u+)  -  Dn  u+)  (p)  =0,  Pe3D, 

or,  rearranging  terms, 


(2.19)  (-a  +  Ka  +  Dn)  u+  =  (K-I)g+, 
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where 

(2.20)  (D  u  )  (p)  :  =  lim_  l—  (Du)  (P)  . 

P-op  p 

Care  must  be  exercised  to  ensure  existence  of  this  limit.  We  deal 
with  this  question  (which  depends  on  both  the  smoothness  of  the 
density  and  the  degree  of  regularity  of  3D)  in  Lemma  4.1. 

As  with  the  integral  equations  derived  in  [ 6  ] ,  we  will  see 
that  the  representation  (2.15)  affords  a  solution  of  the  exterior 
boundary  value  problem  (the  Robin  problem  in  the  present  case) 
provided  u+eL£ (3D)  satisfies  the  pair  of  boundary  integral 
equations.  In  fact,  as  we  shall  see  below,  solutions  of  the  system 
(2.14)  and  (2.19)  will  necessarily  be  more  regular  than  merely 
square-integrable . 

Moreover  we  will  see  that,  as  in  [6  ],  this  pair  of  boundary 
integral  equations  is  redundant  unless  the  homogeneous  adjoint  of 
equation  (2.14)  has  non-trivial  solutions.  In  such  a  case,  how¬ 
ever,  the  pair  of  integral  equations  will  have  a  unique  solution; 
it  is  this  result  we  wish  to  establish  in  the  following  pages. 

Note,  however,  that  unlike  the  cases  discussed  in  [6  ] ,.  the  second 
of  these  equations  is  an  operator  equation  of  the  second  kind. 

The  spectrum  of  K  is  clearly  of  vital  importance  in  discussing 
the  solutions  of  these  equations  and,  since  K  may  be  considered  as 
an  operator-valued  function  of  the  parameter  k,  it  is  convenient 
to  treat  separately  the  operators  (-K)  and  (+K)  rather  than  intro¬ 
duce  another  parameter,  X,  in  order  to  obtain  the  standard  form 
(I-XK) .  To  this  end,  we  adopt  the  following  definition.  Those 
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oo  .  . 

real  values  of  k  for  which  the  equation  w-Kw  =  0  has  nontrivial 

o 

solutions,  w,  will  be  called  characterisitc  values  of  K. 

o  o 

Similarly,  those  real  values  of  k  for  which  w  +  Kw  =  0  has  non- 

o 

trivial  solutions,  w,  will  be  called  characteristic  values  of 
(-K) .  We  note  that  if  k  is  a  characteristic  value  of  K,  then  it 
is  also  a  characteristic  value  of  K*  [[5lChap.  XIII,  §1.3],  and 
trivially  it  is  also  a  characteristic  value  of  K  and  K*.  Similarly, 
(-K) ,  ( — K* ) ,  ( — K> , (-K*)  have  the  same  characteristic  values,  which 

may  differ  from  the  preceding.  Further,  if  w  is  a  nontrivial 

o  o  2  0-2 

solution  of  w  +  Kw  =  0,  then  w  is  a  nontrivial  solution  of  w  +  Kw=0 , 

and  &*+K *&*  =  0  implies  w*+R*w*  =  0.  hereafter,  an 

elevated  index  zero,  for  example  w,  will  indicate  that  the  function 

concerned  is  a  solution  of  a  homogeneous  equation. 

The  values  of  k  for  which  there  exist  nontrivial  solutions 

of  the  homogeneous  Helmholtz  equation  in  D_  with  vanishing  normal 

derivative  on  3d  will  be  called  eigenvalues  of  the  interior 

Neumann  problem.  Similarly,  eigenvalues  of  the  exterior  Neumann 

and  exterior  and  interior  Dirichlet  problems  are  those  values  of 

k  for  which  there  exist  nontrivial  solutions  of  the  corresponding 

boundary  value  problems. 

It  is  well  known  tI-14]  vol.  IV,  §229]  that  there  are  no 
eigenvalues  of  the  exterior  problems ,  but  there  are  eigenvalues 
of  the  interior  problems,  and  furthermore,  these  eigenvalues  are 
real  (Stakgold  [  [15]  vol.  II,  p.  137])..  The  relationship  between 
characteristic  values  of  K  and  eigenvalues  of  the  interior  problems 
is  clarified  in  the  following  theorem,  which  appears  in  this  form 
in  [5 ,  p.  222] . 


Theorem  2.1: 


(A)  k  is  a  characteristic  value  of  K  if  and  only  if  k  is 
an  eigenvalue  of  the  interior  Neumann  problem. 

(B)  k  is  a  characteristic  value  of  (-K)  if  and  only  if  k 
is  an  eigenvalue  of  the  interior  Dirichlet  problem. 

The  interested  reader  is  referred  to  [ 6  ]  for  the  proof  of 
this  result.  We  remark  only  that  an  immediate  consequence  is 
that  if  k  is  an  eigenvalue  of  the  interior  Dirichlet  problem, 
then  k  is  a  characteristic  value  of  (-K)  and  hence  ( — K* ) . 

To  conclude  this  section,  we  wish  to  point  out  that  the 
analysis  presented  here  is  applicable  to  the  scattering  problem 
for  imperfectly  reflecting  surfaces.  Under  appropriate  hypotheses, 
among  which  is  one  requiring  the  radius  of  curvature  to  be  large 
relative  to  the  penetration  depth,  the  exact  transmission  condi¬ 
tions  may  be  approximated  by  a  boundary  condition  of  the  third 
kind  with  purely  imaginary  a  which  may  interpreted  as  an  acoustic 
impedance  [13] .  The  analogous  problem  in  electromagnetics 
is  discussed  by  Senior  [12].  It  is  possible,  then,  to  reduce 
the  scattering  problem  to  an  equivalent  pair  of  boundary  integral 
equations.  These  equations  differ  from  equations  (2.14)  and  (2.19) 
only  in  their  right  hand  members  and,  as  we  shall  see,  for 
sufficiently  smooth  incident  fields  the  analysis  of  Section  4  will 
be  valid  for  the  scattering  problem. 

Specifically,  we  may  proceed  as  follows*  Denoting  the  incident 
and  scattered  fields,  respectively  by  u^  and  us  and  recognizing 
that  the  incident  field  must  be  a  solution  of  the  Helmholtz  equation 
on  the  interior  domain,  we  may  use  the  Helmholtz  representations 
(2.8)  and  (2.9)  to  show  that,  on  the  boundary  3D 
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(2.21)  uS  *  S(3us/3n)  -  K*us 

(2.22)  u1  =  k*ua  -  S(3ui/3n). 

Consequently , 

u  =  us  +  u1  =  2U1  +  S ( 3u/3n)  -  K*u. 

Invoking  the  boundary  condition  3u/3n  +  cu  =  0 ,  this  last  relation 
may  be  rewritten  as 

(2.23)  (I  +  Sa  +  K*)u  =  2uX 

where,  we  emphasize,  u1  is  the  known  incident  field. 

Likewise,  since 

(2.24)  uS(P)  =  (1/2)  [S(3us/3n)  -  DuS]  (P)  ,  PeD+, 

and 

(2.25)  uX(P)  *  (1/2)  [Du1  -  S  (  3ui/3n)  ]  (P)  ,  PeD  _, 

we  have,  taking  normal  derivatives  and  taking  account  of  the 
jump  conditions,  that  for  pe3D, 

3us/3n  =  ( 1/2) [ 3uS/3n  +  K  3uS/3n  ]  -  (1/2)  D_us, 

tr  r  XT 

3ui/3n  =  (1/2)  [  3uV3n  -K  3uX/3n  ]  +  (1/2)D  u1. 
p  p  p  n 


and 


It  is  easy  to  check,  from  these  last  two  relations,  that  since 


3u/3np  =  -a  u  on  3D,  the  total  field  u  must  satisfy  the  equation 

(2.26)  (-a  +  K  o  +  Dju  *  2  (3u1/3n  ). 

n  p 

Equations  (2.23)  and  (2.26)  are  thus  the  particular  forms  of 
(2.14)  and  (2.19)  applicable  to  the  scattering  problem. 

3.  Auxiliary  Results 

In  this  section  we  collect  some  known  results  and  their 

extensions  as  well  as  some  new  results  which  will  be  used  in  the 

subsequent  discussion.  For  more  details,  the  interested  reader 

should  consult  the  cited  references;  we  have  attempted  to  select 

not  original  references  to  these  particular  theorems  but  rather 

those  which  are  readily  accessible. 

We  begin  with  the  result  alluded  to  above  which  we  will 

refer  to,  in  the  sequel,  as  the  regularity  theorem.  The  particular 

version  of  the  results  presented  here  is,  on  the  one  hand,  a 

specialization  to  weakly  singular  surface  integrals  of  a  theorem 

of  Mikhlin  [10]  for  closed,  bounded  sets  in  tn  which,  on  the 

other  hand,  is  a  minor  generalization  of  Mikhlin' s  theorem 

because  of  the  occurrence  of  the  factor  u  which  lies  in  L  (  D) . 

00 

Theorem  3.1: 

Let  3D  be  a  Lyapunov  surface  (or  a  finite  number  of  dis¬ 
joint  Lyapunov  surfaces)  .  If  A;  3Dx3D-M£  and  f:  3D-HE  are  continuous, 
ueL^OD)  and  a<2,  then  any  solution  of 


t 
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(3.1) 


U^P0) 


A<P,P0> 

Ra(P/P0) 


y(p)  u  (p)  dS 

P 


f<p0» 


which  belongs  to  I^OD)  is  also  continuous  on  3D. 

For  a«y  e>0,  choose  n(t)  to  be  a  continuous,  real¬ 
valued,  non-increasing  function  of  t>0,  with  the  property  that 
n(t)»l,  0<t<  e/2,  and  n(t)=0,  t>e.  Then,  setting 


(3.2)  Kx  (p,pQ) : 


and 


A(P>Pq)  n  (R(p,p0)) 
Ra(P,P0)  ~~ 


(3.2)  Kj (p  *Pq)  J 


A(P*P0)  tl-n (R(p ,pQ) ) ] 
Ra(p»P0) 


where  K2(p,p0)  is  continuous  on  3Dx3D,  the  integral  equation  (3.1) 
may  be  written  as 

(3.4)  u(p0)  -  (KlUu) (pQ)  *  F(p0) 


where 


(3.5)  (K^yu) (pg) : 


r 


Kl(pfPg)  y(p) 


u(p) 


and 


(3.6)  F(p-)  -  f(p  )  + 


3D 


K2fp'p0)w(p,u{p*  dst 


The  function  F  is  continous  since  f  is  continuous  and 


(3.7) 


3D 


[K2 (p,pQ) -K2 (p,Px) ] y (p) u(p) dSp  j < 


( 


>  3D 


K2(p,pQ)-K2(p,p1)  r  dS 


1/2 


yu  | 


L2(3D) 


where  K2  is  continuous  in  both  variables.  Note  that  hue  L2(3D) 
since  ue  L2  (3D)  and  ueL,,,  (3D). 

Now  observe  that,  with  the  Schwartz  inequality  and  the  fact 
that  n<l,  and  writing  R  for  R(p,pQ) , 


(3.8) 


I  |KlUu|  | 


2L2(3D)“ 


-  3D 


A(p,P0)n(R)u(P) 

- 7572 - 


2 


f 

r 

3D 

i 

a(p,p0) 1 2n2  (R) |y(p) ,2 


3D 


dS 


UMlfds 

3D  Ra  P 


<  c 


f  i 

t  dS 

!  — E 

f  dS  \ 

'  3D  * 

'  3DoB  (pn)  Ra 

J  3D  R“  PJ 

where  c‘ 


C(3Dx3D) 
o 


L  (3D)  and  B£  (pQ)  is  a  ball  of 


radius  e  centered  at  Pg,  outside  of  which  n(R(p,p.)«o. 
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If  d  is  the  radius  of  the  Lyapunov  sphere  associated  with  3D 
and  e<d,  we  have  the  following  estimates: 


f  dS 

(3.9}  |  — £  <  c' 

J3DnB  (pQ)  ~ 


0l-“dp  ,  £_SL. 

2-a 


.  2-a 


for  some  c'  independent  of  pQ  and  e  and 


(3.10) 


dS 


Ed 


3D  R 


a 


dS. 


Ed  + 


dSPfl_  ,  c'd2“a, 


3DnBd(pQ)  Ra  ^D\3DnBd(pQ)  r' 


,a  —  2-a 


where  S  is  the  surface  area  of  3D. 

Employing  these  estimates  in  (3.8)  we  find  that 


(3.11)  IlK.uu 


2  <  C2C*  +  4  1  !u 

r  (  ryv\\  C  ^  o”  /n, 


l*'*11  L-  ( 3D)  -  w  w  T=a  v2^a  d~  "2 


L~(3D) 


Thus  for  e  sufficiently  small 


(3.12)  IIV'I^OD)  Kl 

and  the  solution  of  (3.4)  in  Lj ( 3D)  is  given  by 

00 

(3.13)  u(pQ)  *  l  ( (K1uJ  nF)  (pQ)  • 

n=0 


Wl  TJ 


16 


In  fact,  this  series  is  uniformly  convergent  hence  u(pQ)  is 
continuous.  This  follows  because  F  is  continuous,  hence  bounded 
on  3D,  and,  using  arguments  similar  to  the  preceeding  (in 
particular  using  (3.9))  we  may  establish 

t  A(p,pn) 

(3.14)  |  (K.yF)  (p0)  H  - — ^n(*)u(p)  F (p)  dS  | 

1  U  J 3D  Ra  p 


<  c  If 


Ra  <  cc 

C£3D>  ^  3DnB£  (pQ )  ~  ^ 


cc '  2 -a , I  , 

t-tt  E  F 


COD) 


Thus  we  have  the  estimate 


(3.15)  | | (KlU)nF| | C(3D)<  e2-“| | (Kju) n-1F| ]c( 


.  fcc 1  2-a.n, i _ | | 

1  (5^T  e  >  Nf|  lc(3D) 


which,  for  s  sufficiently  small,  establishes  the  uniform  conver¬ 
gence  of  (3.13)  and  hence  the  continuity  of  u. 

The  next  three  results  may  be  found  in  the  book  of  Gunter  [ 4  ] 
and  relate  to  the  behavior  of  single  and  double  layer  potentials . 
The  results  on  potentials  with  integrable  densities,  due  originally 
to  Fichera  [3],  are  proven  in  (4]  for  the  case  k*0  but  remain  valid 
for  k^O  since  the  singular  part  of  the  kernels  concerned  remains 
unaltered.  The  first  of  these  results  establishes  continuity 
properties  of  single  and  double  layer  potentials. 
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Theorem  3.2: 

If  3D  is  a  Lyapunov  surface  (or  a  finite  number  of  disjoint 
Lyapunov  surfaces)  and  if  ueL^OD)  then  the  single  layer  potential 
with  density  u  is  H31der  continuous  in  the  entire  space,  while 
the  double  layer  potential  with  the  same  density  is  Holder  con¬ 
tinuous  on  3D.  Moreover,  both  are  infinitely  dif ferentiable  off 
the  boundary. 

Remarks :  The  statement  of  this  theorem  may  be  summarized 
by  the  assertion:  ueL^OD)  implies  (a)  Sue  CQ  a(*3)  nC^d^SD)  and 
(b)  K*u  £  Cn  ( 3D) ,  Dye  C  d3\3D)  .  The  statement  here  differs 
from  that  appearing  in  Gunter  in  that  the  assumption  that  yeLs(3D) 
replaces  the  assumption  that  u  is  bounded  on  3D.  The  estimates 
necessary  for  the  proof  are  identical  in  the  present  case  to  those 
appearing  in  Gunter. 

Turning  our  attention  to  the  jump  relation,  we  may  state  the 
the  following  result.  Recall  that  if  u  is  integrable,  then 
almost  every  point  of  its  domain  is  a  Lebesgue  point. 

Theorem  3.3: 

If  3D  is  a  surface  as  described  above,  ueL1(3D)  and  if  pe3D‘ 
is  a  Lebesgue  point  of  u ,  then 


(3.16) 


.  .  3 

lim  «, 

_  -  3n 

P-p  p 


(Su)  (P) 


(-U+Ku) (p) / 


(3.17) 


(Su)  (P)  -  (U+Ku)  (p)  , 


L8 


and,  for  the  double  layer, 

(3.18)  lim  (Du)(P)  =  (u+K*u)(p), 
P-*-p" 


and 


(3.19)  lim  (Du) (P)  =  (-u+K*u)(p). 
P-p+ 


Moreover,  the  limits  in  each  of  the  four  cases  is  integrable. 

We  complete  our  discussion  of  the  properties  of  single  and 
double  layer  potnetials  with  the  following: 


Theorem  3.4: 


Let  3D  be  a  Lyapunov  surface  (or  a  finite  number  of  disjoint 

A  «v 

Lyapunov  surfaces)  with  index  1  (i.e.  np  •  np  <  ER(p,p2)  where 

y.  1  2 

nPi,  i=l,2,  is  the  exterior  unit  normal  to  3D  at  the  point  p^  and 


E  is  a  constant  independent  of  PlfP2  and  let  u  be  continuous  on 


3D.  If  the  double  layer  potential  (2.4)  possesses  one  of  the 
derivatives 


(3.20)  Du  (p)  :  =  Duu(p)  or  Du  (p)  :  =  Dnu  (p) 

P  P 

then  both  derivatives  exist  and  are  equal  at  the  point  pe9D. 
Their  common  value  will  be  denoted  by  Dnu(p). 
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Next,  we  present  a  version  of  the  divergence  theorem  which  is 
essential  to  what  follows.  This  treatment  allows  the  development 
of  an  L2  theory  of  surface  potentials  which  is  of  interest  for  its 
own  sake  as  well  as  the  present  application.  It  is  closely  allied 
to  the  treatment  of  the  volume  integral  as  an  improper  integral 
presented  in  Smirnov  [14] .  We  remark  that  an  alternative  generalized 
divergence  theorem  involving  potentials  with  continuous  rather  than 
square  integrable  densities  and  based  on  the  interpretation  of  the 
normal  derivative  as  a  distribution  or  "boundary  flow"  is  available 
HI /  18]  and  has  been  used  in  an  integral  equation  approach  to  the 
Neumann  problem  for  the  Helmholtz  equation  [7] .  That  treatment  does 
not  address  the  non-uniqueness  problem  at  characteristic  values  of 
the  integral  operator  in  contrast  to  the  present  approach. 

Theorem  3.S: 

Let  3D  be  a  Lyapunov  surface  of  index  1  and  let  {S^-}  be  a 

family  of  "parallel"  surfaces  in  D_  such  that  (i)  for  fixed  5  , 

S'  is  the  image  of  3D  under  a  continuously  differentiable  trans- 
a 

formation  F.":  3D-*-S  (ii)  P7(p)  :*  P**  P+0(S)  and  if  p(u,v)  is 

n  »P  . 

a  parameterization  of  3D,  Tu  =  fu  +0(6)  and  av'* 

8r  +  0(5);  (iii)  dSp  =  j  dSp  +  0(5.)  for  every  £c3d  where 

£6  =  F'(E). 

Now,  if  $eC^(D_)  is  a  vector-valued  function  with  limiting 
values  in  L2 ( 3D)  by  which  is  meant 

(3.21)  |*(PT(P))-*(P)  i2dS  =0(<5°) 

3D  6  p 

then 

(3.22)  n  •  $  (P)  dS_  -  7-*(P)dx  . 

j3D  P  ^  1 D  p 
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Remark  ;  The  order  symbol  0(6)  denotes  a  vector  whose  magni¬ 
tude  decreases  uniformly  with  respect  to  6.  While  it  is  not 
necessary  to  require  that  iP^-pj  be  independent  of  p  for  fixed  5, 
we  do  require  that  the  surface  not  coincide  with  3D  at  any 
point,  i.e., 

(3.23)  inf  |F*(p)-q|  >96 
p,qe3D 

for  some  9>o,  independent  of  6. 


Proof:  Since  $eC1(D_)  and  D_,  the  divergence  theorem 

in  its  classical  form  yields 


(3.24) 


n  •  $(PJ  dS  = 
P  6 

6  6 


V  •  $  (P)  dx  . 


where  denotes  the  interior  of  S^.  Consider,  first,  the 

integral  on  the  left  hand  side.  The  surface  3D  may  be  decomposed 

N 

into  a  finite  number  of  patches,  3D=  (J  E.  where  each  E.  ,  i=l,...,n, 

i=l  1  1 

lies  in  a  Lyapunov  sphere.  Let  Z^  =  F^(E^).  Then 


(3.25) 


*-4v4,p*,av  Ji 


Points  on  Z^  may  be  represented  in  local  coordinates  (£rn»C(Crn)) 
where  Z.  is  the  image  of  some  domain  D . r  in  the  (£, n) -plane. 
Similarly,  (£,n)  may  be  chosen  to  parameterize  Z^.  An  element  of 

on 


7T~fT 


area  on  Z.  is  dS  =  /1+  c  c  d^dn  whereas  an  element  of  area 
1  P  F  xl 

s  - - — , — 

Z..  is,  using  condition  (ii)  ,  dSp  =[/l+C  + c  +  0(6)]d£dn  or 
16  ’  *  5  E,  n 

dS  =dSp  +0(6). 


Similarly,  on  I.,  n 

^  r 


-?n  1)  -  - 

,  ,  n_  =  - — — — —  ,  while  on  Z.et  =n_+0(6) 
i  p  - = — * -  io  e  .  p 

/ 


Hence,  with  $(P)  and  $ (F^ (p ) ) eL2 ( 3D) ,  we  have 


(3.26)  n  •  $  (P  .)  dS  =  n^  •  *  (PT  (p )  )  dS  +0  (6) 

y  Pr  0  P  JyP  6  P 


Zi«  ' 5 


5  JSi 


Also,  in  light  of  formula  (3.21),  we  have  the  estimate 


(3.27)  |  f  n  *$(F"(p))dS  _  [  n  -$(P)dS  |< 

Jr.  P  6  P  Jr  P  P1- 


■  |*(F"(p))-*(p)  |2dS  \1/2+0(6)  <0  (5)  . 

i  J  0  PJ 


Therefore 


(3.28)  lim  n  •  4>(PJdS  =  n  -<t>(p)dS 

«*o  Jris  P6  s  P,  Jr.  p  P' 


6  i 


Thus  the  limit  on  the  left  hand  side  of  (3.24)  exists  as  S~*-0 , 
is  equal  to  the  left  hand  member  of  the  relation  (3.22)  and  hence 
the  limit  of  the  right  hand  side  in  relation  (3.24)  exists  and 
equals  |3D  np  •  4>  (P)  dSp .  Thus  the  integral  D  V*$(P)dSp  has  a 
meaning  as  an  improper  integral  and,  in  this  sense,  we  have 


established  the  equality  (3.22). 
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Remark ;  The  parallel  surfaces  used  here  are  similar  to, 
but  more  restrictive  than  those  discussed  by  Burago,  Maz'ja  and 
Sapozhnikova  in  [1] .  The  differentiability  of  the  transforma¬ 
tion  ensures  the  convergence  of  normals  and  surface  elements 
and,  while  the  assumption  of  such  convergence  may,  in  fact, 
insure  differentiability  of  the  transformation,  we  have  not 
established  this  assertion.  We  also  wish  to  point  out  that  the 
construction  of  parallel  surfaces  given  by  Smirnov  in  [14]  does 
not  guarantee  that  the  normals  converge  unless  the  surface  3D  is 
in  fact  C2,  i.e.  the  normal  is  differentiable  rather  than  merely 
Lipschitz  continuous. 

Theorem  3.5  has  two  immediate  consequences  which  we  state 

here . 

Corollary  1:  If  we  C2(DJ  and  if  w  and  have  limits  in 

L2 ( 3D)  from  D_  in  the  sense  of  Theorem  3.5,  then 

D  (|Vw(P)  |2+w(P)  a  w(P))dfp. 

Corollary  2:  If  w<e  C2(D+)  ,  if  w  and  ~  have  limits  in  L2(3D) 
from  D+  in  the  sense  of  Theorem  3.5,  and  if  BR  is  a  ball  of 
radius  R  containing  D  in  its  interior,  then 


(3.29) 


w(p) 


(3.  30) 


3w 


3D 


w(p)1H-(p)dS  =  - 


(  |  Vw(P) 


VBR 


•+w(P)  A  w(P)  )  dt  + 
P 


3B 


w(p)!£(p)d 


The  first  of  these  corollaries  follows  immediately  from 
Theorem  3.5  by  identifying  $  in  that  theorem  with  the  function  wVw. 
The  second  follows  from  the  first  by  applying  the  statement  to  the 
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interior  domain  bounded  by  3D  and  3B 

R  • 

We  conclude  this  section  with  a  uniqueness  theorem  for  the 
exterior  Robin  problem  for  the  Helmholtz  equation.  The  proof 
is  based  on  Rellich's  lemma  and  was  given  by  Leis  [9  ]  in  the 
classical  case  for  real,  continuous  a.  The  extension  presented 
here  takes  into  account  complex-valued  a  in  Lm(3D)  and  considers 
a  more  general  class  of  solutions.  Specifically,  we  consider 
solutions  of  the  Helmholtz  equation  in  D+  which,  together  with 
their  normal  derivatives  take  on  boundary  values,  square- 
integrable  on  3D,  in  the  sense  of  Theorem  3.5.  This  notion  is 
precisely  the  notion  of  generalized  boundary  conditions  described 
by  Miranda  [$;  §29].  We  refer  the  reader  to  that  work  for 
a  discussion  of  the  history  of  this  interpretation  of  boundary 
conditions  for  elliptic  equations.  It  is  however  vital  to  the 
results  of  section  IV  to  know  that  single  and  double  layer 
distributions  assume  their  boundary  values  in  this  generalized 
sense  which  we  state  in  the  following: 

Theorem  3.6:  If  3D  is  a  Lyapunov  surface  (or  a  finite  number  of 
disjoint  Lyapunov  surfaces)  and  yeL,  (3D)  then  Sy ,  Dy ,  and 

3n— 

assume  boundary  values  in  L2 (3D)  in  the  sense  of  Theorem  3.5. 

The  proof  of  this  theorem  is  deferred  to  the  appendix.  Through¬ 
out  the  remainder  of  this  paper,  the  boundary  value  problem  will 
be  understood  in  this  generalized  sense  and  the  statement  of 
the  boundary  value  problem  (2 . 10) - (2 . 12)  should  be  so  interpreted. 
With  this  understanding  we  state 
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Theorem  3.  7 : 

Let  3D  be  a  Lyapunov  surface  of  index  1,  let  kel  with 
0<arg  kor  (lmk>0)  and  let  cteL^OD)  with  arg  k<  arg  o<_ arg  k+-rr. 
If  weC2(D+)  is  such  that 

(3.31)  (A  +k2)  w(P)  =0  #  PeD+  ,  - 


3w 


(3.32)  p*  -  ikw  =  o(l/r),  r=|P| 


and 


^  w 

(3.33)  — v(p)  +  aw(p)=  0,  pe3D  , 

3nP 


then  w=0  almost  everywhere  in  D+u3D. 


Remark :  We  repeat  for  emphasis  that  the  boundary  condition 

(3.33)  is  to  be  interpreted  in  the  generalized  sense  and  the 
equality,  consequently,  holds  only  almost  everywhere  on  3D. 

Proof:  Let  BR  be  a  ball  with  boundary  3BR,  center  in  D_  and 
radius  R  sufficiently  large  so  that  D  cB_.  Assuming  that  w 
satisfies  the  hypothesis  of  the  theorem,  we  may  apply  Green’s 
theorem  to  the  function  w  and  its  complex  conjugate  w  in  the 
region  D+nBR  obtaining,  using  Theorem  3.5  and  the  relation  (3.30) 
the  integral  identity 


(3.34) 


(wA 

nB„ 


w-wa  w) df= 


(w- 


3B. 


3w 

In 


-  w 


-  3w 
3n' 


)  dS 


3D 


(”TO 


-  5  !£>ds- 
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The  radiation  condition  (3.32)  shows  that 


(3.35)  (|2  -  ikw)  (~*iiw)=o(I2) 

from  which  we  may  deduce  that,  on  3BR, 


(3.36) 


Substituting  this  result  in  the  identity  (3.34)  and  making  use  of 
the  fact  that  both  w  and  w  are  solutions  of  the  Helmholtz  equation 
and  satisfy  homogeneous  boundary  conditions  (equations  (3.31  and 
(3.33)  and  their  conjugates)  we  obtain 


(3.37) 


7  ~2 

(k‘-k  1 


|w|2dt-  jjjaf  {  {'S  +l*w|2  +  (lmk)  (w|£+  w|£)}< 

D+nBR  3BR 


(a-o)  IwrdS  +  o(l) 


3D 


However,  Theorem  3.5  may  be  applied  to  wA  w+wA  w  and,  again  using 
the  boundary  conditions,  we  obtain 


(3.38) 


(wi£+wi£)  dS=  f  (2|7w|2-(k2+k?|w|2)dT-  f (o+ff)  |w| 
3Br  3n  3n  JD+nBR  i  3D 


2 


dS. 


Substituting  this  result  in  (3.37)  and  making  use  of  the  facts 


(3.39)  i(Rek) (k2-k2) - (Imk) (k2+k2) =2 | k | 2 (Imk) 


and 
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(3.40)  (Rek) (Ima) - (Imk) (Rea) *  Im(ka) 


we  obtain 


(3.41) 


(||£|2+|kw|2)dS+2(Imk) 

3BR 


( |kw|2+| Vw|2)dt 


VBR 


+  2  (  (linker)  |  w|  2dS  -  o(l)  . 

>  3D 

All  the  terms  on  the  left  are  easily  seen  to  be  non-negative  if 
(lmk)>0  and  (lmko)>0,  this  last  condition  being  equivalent  to  the 
condition  arg  k^arg  a<arg  k+7r.  it  follows,  then,  from  Rellich's 
Lenuna  that  w=0 . 

Note  that  the  space  in  which  w  lies  allows  us  to  conclude 
that  w=0  in  D+  whereas  we  can  assert  only  that  w=0  almost  every¬ 
where  on  3D.  It  should  also  be  noted  that  the  classical  solution 
space,  C2 (D+) nC^ (D+u 3D)  is  a  subspace  of  the  space  of  functions 
specified  in  the  statement  of  the  theorem  hence  in  the  classical 
case  w=0  everywhere  in  D+u3D. 

4.  Existence  and  Uniqueness  of  Solutions 

The  formal  calculations  of  section  2  have  given  us  a  pair  of 
boundary  integral  equations  for  a  function  ueL2(3D),  namely  the 
equations 

(4.1)  (I  +  K*  +  Sa)u  =  Sg+ 

(4.2)  (-a  +  Kc  +  Dn)  u  *  Kg+  -  g+. 
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Our  object  in  this  section  is  to  demonstrate  the  existence  of  a 
unique  solution  of  this  pair  of  equations  and  to  establish  the 
equivalence  of  this  problem  to  the  original  boundary  value  problem. 
The  proof  will  proceed  by  means  of  a  series  of  lemmas.  Through¬ 
out  this  section,  unless  explicitly  mentioned  to  the  contrary, 

(e.g.  Lemma  4.3)  the  functions  g+  and  a  will  be  assumed  to  lie 
in  the  space  L^OD)  .  Furthermore,  the  boundary  3D  will  be  assumed 
to  satisfy  the  hypothesis  of  theorem  3.5. 

We  remark,  first,  that  by  applying  the  regularity  theorem 
(Theorem  3.1)  to  any  solution  of  equation  (4.1) ,  we  may  conclude 
that  any  L2-solution  of  equation  (4.1)  must,  in  fact,  be  continuous. 
Note  that  the  regularity  theorem  is  applicable  here  since,  under 
the  hypotheses  on  g+,  the  function  Sg+  is  continuous  everywhere. 
(Theorem  3.2),  and  the  kernels  of  the  integral  operators  R*  and  Sc 
are  products  of  weakly  singular,  essentially  bounded  and  continuous 
functions.  In  fact  it  is  to  ensure  the  continuity  of  Sg+  and  the 
subsequent  applicability  of  the  regularity  theorem  that  we  restrict 
g+  to  lie  in  rather  than  L2 .  As  we  shall  see  below  the  fact 
that  1*2  solutions  of  (4.1)  must  also  be  continuous  will  be  of 
crucial  importance  in  ensuring  that  certain  functions  used  in  the 
construction  of  the  required  solution  of  the  system  (4.1) -(4.2) 
will  lie  in  the  domain  of  the  operator  Dn.  Indeed,  we  will  begin 
with  the  following  basic  result  relating  the  set  of  solutions  of 
equation  (4.1)  and  the  domain  of  this  operator. 

Lemma  4,1:  Suppose  that  the  wave  number  k  is  not  an  eigen¬ 
value  of  the  interior  Dirichlet  problem  for  the  Helmholtz 
equation.  Then,  for  any  solution  u  of  equation  (4.1)  in  L^OD) 

D  u,  and  hence  D*u,  exists  and,  moreover,  D-u  ■  D^us-D  u . 

11  n  n  n 
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Proof:  Let  u  be  any  solution  of  the  equation  (4.1) .  In 
light  of  the  theorem  of  Gunter  (Theorem  3.4  above)  it  will  be 
sufficient  to  show  that  D~u  exists,  at  least  almost  everywhere 
on  3D. 

To  this  end,  define  the  function  v  on  D_  by  the  equation 

(4.3)  v(P)  :=  (Du  +  Sau) (P)  -  Sg+(P),  PeD_  . 

Then  v  takes  on  boundary  values  in  the  sense  of  Theorem  3.5. 

Since  the  function  u  is  a  solution  of  the  integral  equation 

(4.1)  ,  the  regularity  theorem,  3.1,  assures  that  u  is  continuous 
on  3D.  Therefore  u  is  bounded  on  3D  and  hence  the  function  au 
is  essentially  bounded,  measurable,  and  so  integrable  on  3D. 
Likewise,  the  function  cru  -  g+  is  essentially  bounded  and 
integrable  on  3D,  and  so  the  single  layer  with  density  au  -  g+ 
is  continuous  everywhere  by  Theorem  3.2.  Recalling  the  form  of 
the  jump  conditions  for  double  layer  with  continuous  density 
(theorem  3.5),  we  have  for  all  pc3D, 

(4.4)  lim  v(P)  *  (u  +  K*u  +  Sa  u) (p)  -  (Sg  )  (p) . 

P-p"  + 

Since,  by  assumption,  the  function  u  is  a  solution  of  the  equation 

(4.1)  we  may  conclude  that  lim  v(P)  =0  for  all  points  pc  3D. 

P-*p" 

The  function  v,  therefore,  is  a  solution  of  the  interior  homogeneous 
Dirichlet  problem  with  zero  boundary  data.  Since  k  is  not  an 
eigenvalue  of  this  problem,  the  function  v  must  vanish  identically 
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in  6  and  so  3v/3n”  =  0  on  3D. 

P 

We  now  rewrite  the  relation  defining  v  in  D_  as  follows: 

(4.5)  (Du) (P)  *  v(P)  +  Stg.  -  au] (P) ,  PeD 

*  ™  • 

Again,  since  density  of  the  single  layer  operator  in  this  last 
relation  is  integrable ,  we  may  use  Theorem  3.3  to  conclude  that 

(4.6)  lim  *  Stg.  -a  u]  (P)  =  [-g,  +  Kg,  +au  -Kau]  (p) 

P-p"  3np 

almost  everywhere  on  3D.  Since  we  have  just  shown  that 
3v/3np  ■  0  for  all  pc 3D,  we  may  conclude  that  Dftu  exists  for 
almost  all  pe3D  which  is  what  we  wished  to  demonstrate.  Moreover 
at  such  points  p, 

(4.7)  D~u (p)  *  (-g+  +  Kg+  +  au  -  Kau)  (p) , 

which  relates  the  solutions  of  equation  (4.1)  to  those  of  equation 
(4.2)  .  More  precisely,  we  have  established  the  following  result. 

Lemma  4.2:  If  k  is  not  an  eigenvalue  of  the  interior  homogeneous 
Dirichlet  problem  and  if  the  function  u  is  a  solution  of  equation 
(4.1)  in  L2 (3D) ,  then  the  function  u  is  likewise  a  solution  of 
equation  (4.2) . 

We  may  now  turn  to  a  discussion  of  the  existence  and  unique¬ 
ness  of  solutions. 


Lemma  4.3:  For  all  choices  of  g+eL2(3D),  the  equation  (4.1) 
has  a  solution  in  L2(3D) . 

Proof:  The  integral  operators  in  (4.1),  K*  and  So,  are 
completely  continuous  on  L2(3D)  (Mikhlin  [10])  hence  Fredholm's 
Alternative  applies  to  (4.1) .  Observe  that  the  operator  So  should 
be  viewed  as  a  composition  of  multiplication  by  the  essentially 
bounded  function  c  which  maps  T>2  ( 3D)  into  itself  and  the  single 
layer  which  is  weakly  singular.  Fredholm's  Alternative  implies 
that  for  any  g+  either  equation  (4.1)  has  a  solution  or  the  homo¬ 
geneous  problem  has  non-trivial  solutions,  in  which  case,  equation 
(4.1)  has  a  solution  if  and  only  if  the  function  Sg+  is  orthogonal 
to  all  solutions  of  the  homogeneous  adjoint  equation  which  may 
be  written  (I+K+crS*)w  *  Oor,  equivalantly ,  (I+K+a§)w  »  0.  In  fact, 
we  will  show  that,  for  all  such  w,  (Sg+,w)  «  0  for  any  choice  of 
g+eL2(3D) . 

To  this  end  note  that  (Sg+,w)  *  (g+,Sw)  and  define  a  function 
v  in  D+  by 

(4.8)  v(P)  :«  (Sw)  (P) ,  PeD+. 

Then,  since  &eCq( 3D)  ,  the  jump  relations  yield,  for  almost  all 
pe3D, 

(4.9)  3v/3n*  (p)  *  (w  +  Kw)  (p) . 

But  the  function  w  satisfies  the  equation  (l+K+aS)$  *  0  and  so 
we  may  write 
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(4.10)  3v/3n*  (p)  =  -oS$  (p)  *  -av(p) ,  for  almost  all  pe3D. 

Hence,  the  function  v  is  a  solution  of  the  exterior  homogeneous 
Robin  problem  which  assumes  boundary  values  in  the  sense  of 
Theorem  3.5  and  so  vanishes  identically  according  to  the  unique¬ 
ness  theorem  discussed  above.  We  may  conclude  Sw  =  0  and  this 
establishes  the  lemma. 

This  last  result,  together  with  the  preceeding  lemma,  shows 
that  if  k  is  not  an  eigenvalue  of  the  interior  homogeneous 
Dirichlet  problem,  then  the  pair  of  integral  equations  (4.1)  - 
(4.2)  has  a  simultaneous  solution.  Moreover,  there  can  be  at 
most  one  such  solution  as  we  now  show. 

Lemma  4.4:  It  UP  L2  0  D)  is  a  solution  of  the  homogeneous 
system  of  boundary  integral  equations 

(4.11)  (I+K*+So)u  »  0 

(4.12)  (-a+Ka+D  ) u  =  0. 

n 

then  u  »  0,  pe3D 

Note  that  since  the  function  u  is  assumed  to  satisfv 
the  equation  (4.11),  we  may  again  invoke  the  regularity  theorem 
to  establish  the  continuity  of  u.  Indeed,  it  may  well  be  smoother 
since  we  assume  as  an  hypothesis  that  it  lies  in  the  domain  of 
the  operator  Dn< 

Proof:  Again,  define  the  function  v  on  D_  by  the  relation 

v(P)  :«  (D  +  Sc)u(P) ,  PeD 
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where  u  is  a  solution  of  the  homogeneous  system  (4 . 11) - (4 .12) . 
Then,  as  before,  v  takes  on  boundary  values  in  the  sense  of 
Theorem  3.5  and,  since  u  is  an  I^-solution  of  equation  (4.11) 
which  has  a  continuous  right  hand  member,  it  is  also  continuous. 
Taking  limits  as  the  point  PeD_  approaches  the  boundary  of  D, 
and  using  the  jump  conditions  for  the  double  layer  operator  with 
continuous  density,  we  have 

(4.13)  v(p)  =  u(p)  +  K*u (p)  +  (Sa)u(p) 

and  so  v(p)  =0  for  all  pe3D  since  u  is  a  solution  of  equation 
(4.11).  Hence  the  function  v  is  a  solution  of  the  homogeneous 
Dirichlet  problem  in  D_.  There  are,  then,  two  possibilities: 
either  (i)  v  vanishes  identically  in  D_,  or  (ii)  v  is  a  non¬ 
trivial  eigenfunction  of  the  interior  homogeneous  Dirichlet 
porblem  for  the  Helmholtz  equation.  In  case  (ii)  the  function 
v  has  regular  normal  derivative  in  the  sense  of  Lyapunov  (see 
Smirnov  [14]  p.  586  and  p.  675) . 

In  the  first  case,  we  must  have  that  (Du  +  Sou) (P)  *  0 
for  all  PeD_.  Now  define  the  function  v+  in  D+  by 

(4.14)  v+(P)  (Du  +  Sou) (P) ,  PeD+, 

which,  like  v_,  takes  boundary  values  in  the  sense  of  Theorem  3.5. 
Using  the  jump  relations  for  the  double  layer  operator  as  the 
point  P  approaches  3D  from  D+,  we  may  write 
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(4.15)  v+  (p)  =  -u(p)  +  K*u  (p)  +  S(au)(p),  pe3D. 

Since,  for  merely  bounded  and  measurable  densities,  the  jump 
relations  for  the  normal  derivative  of  the  single  layer  operator 
hold  only  almost  everywhere  on  3D,  we  may  argue  in  a  manner 
similar  to  the  proof  of  Lemma  4.1  and  conclude  that,  for  almost 
all  pe3D, 

(4.16)  3v+/3n*  =  (au) (p)  +  (Kau) (p)  +  Dnu(p) . 

Note  that  since  k  is  not  an  eigenvalue  of  the  interior  Dirichlet 
problem.  Lemma  4.1  assures  the  existence  of  Dnu(p) .  But,  from 
equation  (4.11) ,  we  have  K*u+Sau  =  -u  and  so,  with  (4.15) , 
v+(p)  ■  -2u(p) ,  pe3D,  while  from  equation  (4.12)  we  have  that, 
with  (4.16),  3v+/3n*  (p)  *  2(au)(p)  almost  everywhere  on  3D. 
Comparing  these  last  two  results,  we  see  that  the  relation 

(4.17)  3v+/3n*(p)  +  (ov+) (p)  =  0 

holds  almost  everywhere  on  3D  and  we  may  invoke  the  uniqueness 
theorem  to  conclude  that  v+  vanishes  identically  in  D+  and 
that  v+(p)  *  0  almost  everywhere  on  3D.  It  follows  since 
v+  *  -2u  that  u(p)  =0  almost  everywhere  on  3D.  But  the  function 
u  is  continuous  on  3D  and  consequently  vanishes  on  3D,  which 
establishes  the  result  in  case  (i) . 
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If  the  second  possibility  obtains,  then  the  function  v  is 
not  identically  zero.  But,  as  before, 

(4.18)  3v/3n*  (p)  =  (Dnu  -  au  +  Kau)  (p) 

for  almost  all  pe3D,  and  so,  since  the  function  u  satisfies 

equation  (4.12),  3v/3n  (p)  =  0  almost  everywhere  on  3D.  But 

P 

then  the  function  v  is  a  solution  to  the  interior  Helmholtz 
problem  which  vanishes  on  the  boundary  and  whose  normal  deriva¬ 
tive  vanishes  almost  everywhere  on  3D.  This  leads  to  a  contra¬ 
diction  and  completes  the  proof. 

It  remains  for  us  to  establish  that,  if  k  is  an  eigenvalue 
of  the  interior  homogeneous  Dirichlet  problem,  then  the  system 
(4.1) -(4.2)  of  boundary  integral  equations  admits  a  simultaneous 
solution. 

We  recall  the  following  result  from  [6,  p.  224  and  p.  277]. 

Lemma  4.5:  There  exists  a  non-trivial  function  w  such  that 

(4.19)  w  +  Kw  =  0 

if  and  only  if  (S$) (p)  =0,  pe3D.  Moreover  such  a  function  exists 
if  and  only  if  k  is  an  eigenvalue  of  the  interior  Dirichlet  problem. 
Finally,  v  is  an  eigenfunction  if  and  only  if  it  may  be  represented 
in  the  form  v (P) *-1/2 (Sw) (P) ,  PeD.  where  $+K&=0 . 

Clearly  such  a  function  $  satisfies  the  boundary  integral 
equation 
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(4.20)  &  +  Kw  +  aS&  =  0, 

and  having  made  this  observation,  we  may  prove  the  following  result 
which  is  analogous  to  Theorem  4.3  of  [g  1. 

Lemma  4.6:  The  function  v  is  an  eigenfunction  of  the  interior 
Dirichlet  problem  if  and  only  if  $  can  be  represented  in  the  form 

o  o 

(4.21)  v  =  Dw*  +  Saw*,  PeD_ 

o 

where  w*  satisfies  the  homogeneous  boundary  integral  equation 

o 

(4.22)  (I  +  K*  +  Sa)  w*  =  0. 

o 

Proof:  Assume,  first,  that  w*  is  a  non-trivial  solution  of 
the  given  boundary  integral  equation  and  hence  with  the  regularity 
Theorem  3.1,  is  continuous.  Define  a  function  v  in  D_  by 

o  o 

(4.23)  v(P)  :=  (Dw*  +  Sow*) (P) ,  PeD_, 

which  takes  boundary  values  in  the  sense  of  Theorem  3.5.  Then, 
as  in  the  proof  of  Lemma  4.1,  the  function  v  vanishes  on  3D  and 
hence  either  vanishes  identically  in  D_  or  is  a  non-trivial  eigen¬ 
function  of  the  interior  Dirichlet  problem.  To  see  that  v  is 
not  identically  zero,  notice  that  if  this  were  the  case  then  its 
normal  derivative  would  vanish  on  3D.  But,  then,  for  almost  all 
pe3D,  we  would  have,  using  the  jump  conditions  as  in  previous 
proofs , 
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_0  0  0 

(4.24)  8v/8rip  (p)  =  (D~w*+Kaw*  -  aw*)  (p)  =  0 


where  D~w*  exists  since  — —  and  —  Saw*  exist  and  with 
n  3n“  dn~o  o 

Theorem  3.4,  which  applies^ because  w*  is  continuous,  Dnw* 

o 

exists.  Thus  the  function  w  is  a  solution  of  the  homogeneous 
boundary  intergral  equation 


o 

(4.25)  (D  +  Ka-a)  w*  =  0 
n 


Now  define  a  function  &  on  D+  by 


o  o 

(4.26)  u(P)  :=  (Dw*  +  Saw*),  PeD+, 


which,  again,  takes  on  boundary  values  in  the  sense  of  Theorem  3.5. 
Again,  using  the  jump  conditions,  we  have  for  P-*-p+, 


(4.27)  u(p)  :=  ( — w*  +  K*w*  +  Saw*) (p) ,  pe9D, 


and 


(4.28)  3S/3n^  (p) 


0  0  0 

(aw*  +  Kaw*  +  Dnw*)  (p) ,  a.e.  on  3D. 


o  o 

But  since  (I+K*Sa)w*=0  we  have  with  (4.27)  that  S(p)  =  -2w*(p)  for 

all  pe9D,  while  from  (4.25)  and  (4.28)  we  conclude  that  3&/an+(p)  = 
o  P 

2(aw*)(p)  almost  everywhere  on  3D.  The  function  S  is  thus  a 

solution  to  the  exterior  homogeneous  Robin  problem  and,  by  the 

o 

uniqueness  theorem,  must  vanish  everywhere.  Thus  w*(p)  =  0 

o 

for  all  pe3D  which  contradicts  the  choice  of  w*.  Hence  v  must 
be  a  non-trivial  eigenfunction. 


fill-  illil  iii'ifiinii 
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Conversely,  assume  that  v  is  an  eigenfunction  of  the  interior 

o 

Dirichlet  problem  and  let  v^,  i=l,2, . . .  ,n,  be  a  basis  of  eigen¬ 
functions.  Thus  there  exists  scalars,  a,,..., a  ,  such  that 
_  in 

o  n  q 

v  a^v^.  lenima  4.5,  to  each  v^  there  corresponds  a  function 

0  o 

such  that  satisfies 

(4.29)  (I+K+aS) wi (p)  =  0  pe3D 

(4.30)  v. (P)  =i(Sw.)(P),  PeD  . 

1^1  — 

4» 

Let  W=span{wj_) .  Note  that  the  i=l,2,.;.,n,  are  linearly  independent 
since  the  v^(p),  i=l , 2 , . . . ,n,  are  chosen  as  a  basis  for  the 
eigenspace.  Hence  dim<$n>=dimW=n  and  Wcker ( I+K+aS) .  But, 

(4.31)  dim  ker  (I+K+aS)  *  dim  ker  (I+K+aS)* 

=  dim  ker  (I+K*+Sa) 

=  dim  ker  (I+K*+So) . 

o 

Let  {w*},  i=l,2,...,rv  be  a  basis  for  ker  (I+K*-rSa)  and  define  for 
each  i=l, . . . ,n 

o 

(4.32)  a.(P)  =  (D+So ) w*i (P) ,  PeD_ . 

2  -  o 

Since  w*  are  non-trivial  solutions  of  ( I+K*+S0)w*=o ,  the  first 

part  of  the  proof  shows  that  the  functions  u^,  i»l...,n  are 
eigenfunctions  of  the  interior  Dirichlet  problem.  Moreover, 
they  are  linearly  independent  for  if  3^,  i=l,...,n,  are  scalars 
such  that  Z6iui=0,  then 
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o 

(4.33)  (D+So) (IBiwt)  =  0 

which  implies,  as  in  the  first  part  of  the  proof,  that 

o 

(4.34)  =  0 

and  hence  8^=0  for  all  i=l,...,n. 

Thus  every  eigenfunction  &  of  the  interior  Dirichlet  problem 
can  be  written  in  the  form 

n  o  n  o  0 

(4.35)  $  =  I  ct  •  <D-*-Sa)  w?  =  (D+Sc)  Z  ct-w*  -  (D+Sa)w* 

i=l  1  1  i=l  1  1 

for  a^,  i=l, ... ,n,  appropriately  chosen  scalars,  which  completes 
the  proof. 

Finally,  we  may  establish  the  following  existence  theorem 
for  solutions  of  the  system  of  boundary  integral  equations 

(4.1) — (4.2) . 

Theorem  4. 1 

Let  aeL^OD)  .  Then  there  exists  a  unique  solution  of  the 
system  of  boundary  integral  equations 

(4.1)  (I  +  K*  +  Sc) u  ■  Sg+ 

(4.2)  (-c  +  Kc  +  Dn)u  =  Kg+  -  g+ 

for  each  choice  of  g.eL  (3D)  . 

+  00 
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Proof:  If  k  is  not  an  eigenvalue  of  the  interior  Dirichlet 
problem  for  the  Helmholtz  equation,  then  Lemmas  4.3,  4.2,  and 
4.4  ensure  that  a  unique  solution  of  the  system  exists  since 
g+eLao(3D)  implies  that  g+eL2  (  3D)  . 

If  k  is  an  eigenvalue,  then  while  Lemma  4.3  still  guarantees 
the  existence  of  at  least  one  solution  of  equation  (4.1),  we 
cannot  show  immediately  that  such  a  solution  will  satisfy  the 
second  equation.  Nor  can  we  make  any  statement  concerning 
uniqueness  since  solutions  to  equation  (4.1)  in  this  case  are 
not  unique . 

Suppose,  then,  that  w^  is  a  solution  of  equation  (4.1) 
and  define  a  function  v  on  D_  by 

(4.36)  v(P)  (Dv^  +  Sov^  -  Sg+)  (P)  ,  PeD_. 

Then,  as  in  the  proof  of  Lemma  4.1,  v  is  either  identically  zero 
or  is  a  non-trivial  solution  of  the  interior  homogeneous  Dirichlet 

problem.  In  either  case,  according  to  Lemma  4.6,  we  may  represent 

o  o 

the  function  v  in  the  form  v  =  (D+Sa)w*  for  some  solution  w*  of 

o 

the  homogeneous  boundary  integral  equation  (I+K*+Scr)  w*=0 .  Hence, 
for  PeD_, 

0  V 

(4.37)  (D+Sa) w* (P)  =  (Dwx)  (P)  +  (Sa)wL(P)  -  (Sg+)  (P)  . 

Now,  define  a  function  w  on  3D  by 

^  2  * 

(4.38)  w  :=  w.  -w* .  _ 
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Then 


(4.39)  w+K*w+Saw  *  w^-w*+K*w^-K*w*+Saw^-Saw* 

%  «•  .a,  f\. 

=  w1+K*w1+Sow1 
=  S*+ 

and  so  the  function  w  satisfies  the  equation  (4.1) .  Moreover, 
for  all  PeD_ ,  with  (4.37), 

o  o 

(4.40)  (Dw  +  Saw) (P)  =  (Dw^  -  Dw*  +  Sow^  -  Saw*) (P) 

*  Sg+(P) 

or  equivalantly , 

(4.41)  (Dw) (P)  «  S(g+-aw)(P),  for  all  PeP_. 

^  2 

The  function  w^  and  w*  are  not  only  -solutions  of  the  non- 
homogeneous  and  homogeneous  integral  equations,  respectively, 
but  are  also,  according  to  the  regularity  theorem,  continuous 
functions  on  9D.  Hence  they  are  both  bounded  there  and,  conse¬ 
quently,  the  functions  w  and  g+  -aw  are  essentially  bounded  and 
integrable  on  9D.  So  again,  we  may  use  Theorem  3.4  to  conclude 


that 
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(4.42)  lim_  3/3n  S(g+  -aw) (P)  =  (-g+  +  Kg  +  aw  -  Kaw) (p) 

P-*.p  p 

almost  everywhere  on  3D.  And  so  we  have  that  D~w(p)  exists 
almost  everywhere  on  3D  and  (Drw  -  ow  +  Kaw) (p)  =  (-g+  +  Kg+) (p) , 
almost  everywhere  on  3D.  Consequently,  the  function  weL2(3D) 
satisfies  the  system  of  boundary  integral  equations.  To  finish 
the  proof,  we  appeal  to  Lemma  4.4  for  the  uniqueness  statement. 

Next  we  establish  the  equivalence  of  the  Robin  Problem  and 
its  boundary  integral  equation  formulation. 

Theorem  4.2: 

Let  g+  and  a  be  in  La(3D) ,lm  k >0 ,  and  Im  ka>0 .  Then  u+  is 
a  solution  of  the  Exterior  Robin  Problem, 

1.  u+  eC2(D+);  u+,  +  eL2(3D)  in  the  sense  of  theorem  3.5 

2.  (A  +  k2)  u+  =  0,  PeD+ 

3u 

3.  yj - iku+  «  o(Vr) 

3u+ 

4.  +  +  au+=g+  almost  everywhere  on  3D, 
if  and  only  if 

(4.43)  u+  «  j  S{g+-au)  -  j  Du,  PeD+ 


where 

(4.1)  (I  +  K*  +  Sa)  u  =  Sg+ 

(4.2)  (-a  +  Ka+  Dn)  u  =  Kg+-g+. 
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Proof:  First  assume  u+  is  a  solution  of  the  Robin  problem 
1-4.  Then  Theorem  3.5  may  be  employed  to  obtain  the  Representa¬ 
tion  Theorem  (2.13)  which  is  the  same  as  (4.43)  above  with  u+ 
replacing  u.  Then  the  formal  calculations  by  which  (2.14)  and 
(2.19)  which  are  identical  with  (4.1)  and  (4.2)  were  obtained  may  be 
repeated.  Note  that  (4 . 1)  follows  from  the  continuity  of  Sg+  and 
Scru+  (Theorem  3.2)  and  the  jump  relation  for  Du+  which  holds 
for  u+eL2OD)  (Theorem  3.3).  Then  the  regularity  Theorem  3.1 

ensures  that  u+eCQ(3D) .  Hence  with  Theorem  3.4  and  the  fact 
3u+ 

that  u+  and  S  have  normal  derivatives  dn  L2 (3D)  the  represen¬ 
tation  (2.13)  may  be  differentiated  and  D  u,  must  exist.  Thus 

n  + 

u+  is  a  simultaneous  solution  of  (4 . 1) and  (4 . 2) hence  by  Theorem  4.1 
is  the  unique  solution. 

Conversely  assume  u  is  the  unique  solution  of  (4.3)  and  (4.2) 

and  u+  is  defined  by  (4.43)..  Clearly  u+  satifies  the  Helmholtz 

equation  and  the  radiation  condition  in  D+. 

The  regularity  Theorem  3.1  ensures  that  u  is  continuous  on 

3D,  hence  Sau  is  continuous  and  has  a  normal  derivative  in  L2  (3D) 

(Theorems  3.2,  3.3,  and  3.6)  .  Moreover  u  must  be  in  the  domain  of 

3u 

+  + 

Dn .  Hence  u+  and  ^  are  in  L2(3D).  Then  taking  limits  as  P-*-p 

we  find  that  on  gp 

(4.44)  u.  =  i  S (g ,  -ou)  +  K*u 

+  2  +  2  2 

and  with  (4.1) 

(4.45)  u+  =  u  on  3D. 


Also 
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(4*46)  3iT  =  k  g+  +  7  Kg+  '  ^  "  I  Kou  "  7  Dnu 
and  with  (4.2) 

3u+ 

(4.47)  3  g+  -  cu  on  3D 

hence 

3u+ 

(4.48)  +  au+  =  g+  on  3D. 

and  u+  satisfies  the  Robin  problem. 

In  conclusion,  we  remind  the  reader  of  the  closing  remarks 
of  Section  2.  There  we  derived,  formally,  boundary  integral 
equations,  of  the  same  form  as  equations  (4.1)  and  (4.2)  ,  but 
with  different  right  hand  members  which  were  suitable  for  describing 
the  total  field  for  the  scattering  problem.  In  that  derivation, 
the  right  hand  sides  represent,  up  to  a  constant  factor,  the 
incident  field  and  its  normal  derivative  respectively.  As  the 
results  of  this  section  do  not  depend  on  the  particular  form  of 
the  function  appearing  on  the  right  hand  side,  we  see  that  for 
sufficiently  smooth  incident  fields  the  analysis  of  this  section 
is  valid  and  we  may  state  the  following  theorem. 

Theorem  4,3: 

Let  u1  denote  an  incident  field  and  u=u1+us  the  total  field 
when  us  denoted  the  field  scattered  by  the  obstacle  D_.  Then,  for 
any  aeL^OD)  ,  there  exists  a  unique  solution  of  the  system  of 
boundary  integral  equations 
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(4.49)  (I+Sa+K*) u  =  2U1 

an1 

(4.50)  (-a+Ka+Dju=2(  /3n  ). 

n  p 

The  equivalence  of  solutions  of  this  pair  of  boundary 
integral  equations  and  the  scattering  problem  follows,  mutatis 
mutandis,  from  Theorem  4.2. 


1 
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APPENDIX  -  On  the  generalized  boundary  values 
of  single  and  double  layers. 

In  Theorem  3.5  the  notion  of  a  family  of  parallel  surfaces 

{ST}  contained  in  D  was  introduced.  We  may  define  another 
o  — 

family  {S*}  contained  in  D+  which  converges  to  3D  in  the  same 
sense  as  {S~}  .  With  the  understanding  that  the  layers  may 

O 

approach  their  boundary  values  from  D+  or  D  we  state  and  prove 
Theorem  3.6. 


Theorem: 

If  3D  is  a  Lyapunov  surface  (or  a  finite  number  of  disjoint 

g 

Lyapunov  surfaces)  and  yeL2(3D)  then  Sy ,  Dy  and  Sy  assume 
boundary  values  in  the  sense  of  Theorem  3.5,  namely; 

A-l  \  |  (Sy)  (F4(p))-(Sy)  (p)  |  2dSn  «o  (6°)  , 

J  3D  *  P 


A-2 


and 


3D 


(Dy)  (P*(p)  )-(Dy)  (p)  |  2dS  =o  (6°)  , 
6  P 


A-S  j  |  (Su)(F6(p))  -  (Sy)  (p)  |  dSp  =o  (6W) 

3  p  P 


on 


Proof: 

Considering  equation  A-l  first 

♦ 

of  a  point  p  on  3D,  i.e.,  P^* 


9 


we  denote  by  Pg  the  image 
(p) .  Recall  that 


A-4  |P5-p|  ;*  R(P5,p)  ■  0(5)  uniformly  in  p. 


46 


With  (2.3)  the  integrand  in  A-l  becomes 


A-5  | 

(Su) (P5)-(Sy) (p) |2=| 

u(q)  [Y(P,r,q) -y(p»q)  Ids  |2 

3D  0  q 

=  _1 

If  ||  ^  />»\ 

r  ikR(P5,q)  ikR(p,q)- 

e  e 

1  dS  |  2 

4ir 

T  1  U iq; 

>  3D 

-  R(P5,q)  R(p,q)  -1 

4 

It  may 

be  shown 

that  for  Xmk>0, 

A-6  | 

ikR(P6,q) 

e 

ikR  (p  ,q) 
e 

(l+|k|D)R(P5,p) 

1  < 

R(P6/q) 

R(p#q) 

R(P5  ,q)R(p,q) 

where  D  is  the  diameter  of  3D.  From  the  definition  of  it 
follows  that  there  is  a  constant  c  such  that, 


A-7  R(P6,p)  <_  c6  . 


Thus  there  is  a  constant  c ,  independent  of  6 ,  such  that  A-5 
gives  rise  to 


A-8 


(  (Su) (P6)-(Su) (p)  i 2  <  c62 


'(  h lM  ldSq _ ' 

/  3D  R(P5 ,q) R(p,q) , 


2 


<  c5 


2 


dS 

q 

R(P5 ,q) R(p,q) 


kl2Lll  dS 

R(P6,q)R(p,q)^ 


47 


* 


where  Schwartz' inequality  has  been  used.  It  also  may  be 
established  that  for  R(P(j,p)<d/2  where  d  is  sufficiently 
small  (depending  on  the  Lyapunov  constants  of  the  surface) 

A-9  5(Pfi,q)  <  K(p,q)  '  Vp'qe3D 

and 

A~10  R(P6,q)  -  ff(Pfi,p)  '  Vp,q£3D  * 

Employing  these  relations  selectively  we  find,  for  any 

a  (0,1)  , 

A-ll  | (SU) (P6)-(Sp) (p) I2  <  c6a 

f  S 

The  term  J  3D  is  an  integral  with  a  weak  singularity 

R  (p,q) 

and  differentiable  density  (=1)  hence  is  a  continuous  function  of 

») 

p  on  3D  while  the  term  [  ^  1 —  dS  is  a  single  layer  with 

}  3D  R(p,q)  q 

integrable  density  hence  is  integrable  (f4],  p.  110)  and  the 
validity  of  A-l  follows. 

Next  consider  A- 2  which  involves  somewhat  different  estimates. 


3D 


-T3- 

1+a 


(p/q) 


I  v  (q)  i 
R(p,q) 


dS 


>  3D 


With  (2.4),  (2.5),  (3.18)  and  (3.19)  we  have 
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A-12  |  (Du)  (Pg) -(D±u)(p)  j2  = 


3  o  * 

W(q)'^'  fY(P,,q)  -Y(p,q)  ]dS^  +u(p)|  ,  P.eD' 
3D  q  0  *1  ~  0 


Introducing  Gauss '  integral  which  has  the  form 


A- 13 


3  3n“  Y  o(P6'q)dSc 

3D  q 


0,  P6eD+ 
•f  1  ,  P^eSD 

2  •  P6eD- 


where 


A-14  YQ(P6,q)  -  -  ,q)  ' 


equation  A-12  may  be  written 


A-15  |  (Dy)  (P^) -(D"u)  (p)  j  2  = 


{u(q)'5rr'[  Y<p*»q)  -Y(p»q)  J 
3D  q  0 


"w(p)lir  ^0(P6,q)-Y0(p,q)]>dS J2 
q 


<  2\  |  u(q)|ir-[Y(P6,q)-Y(p,q)-Y0(P5,q)+Y0(p,q)  ]dS  |2 

<3  H 


3D 


+  2 


3d  Iu(q)-M(p)  j|jj-  [YQ(P5,q)-Y0(p,q)]dSq|2. 


The  triangle  inequality  was  employed  in  obtaining  the  last 
inequality. 
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til 


It  may  be  shown  that  for  lmk>0 


A~16  l|n“[Y  (p6'9)  -Y0(P,s*q)+Y0(p,q)  -YCp*q)  1 1  = 
q 


=1_  3 _ 

2*  3nq 


ikR(P.  ,q)  ikR(p,q) 

e _ 2±.  _  e _ -1 

R(P5,q)  R(p,q) 


< 


c  6 _ 

R(P^ ,q) R (p ,q) 


where  c  is  a  constant  which  depends  on  k  and  the  diameter  of  3d 
but  is  independent  of  p  and  5.  Hence  the  first  term  on  the  right 
hand  side  of  A- 15  may  be  estimated  as 


A-17  2 


u(q)‘3rT[Y(P6'q)  -Yo(P6'q)+Yo(p'q)  “Y(p'*>  IdSl2 
3D  q  M 


1  C,  <52 

-  U  (q-}-l  ds 

-L 

.  . 

3D  R{P  ,q)R(p,q)q 

which  is  exactly  the  same  as  the  estimate  A-7  hence  the  subsequent 
analysis  applies  here  as  well.  The  second  term  on  the  right  of 
A-15  is  analysed  with  the  help  of  the  following  estimate  which 
may  be  shown  to  hold  on  a  Lyapunov  of  order  a: 


A-18 


,9Wq>  9Yo(p'q> 

‘■SrT  "  IxT 
q  q 


c1r(p(S,p)  c2R(P5,P) 

ir(Pg7q)  R(p6,q)  R2(p,q) 


where  c^  and  c2  are  constants.  Then 
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A-i9  lj3D[u(q)-U(p)]  tY0(P5,q)-Y0(p,q)]  dSq  | 


<_  (S  [(  I  U  (q)  -u  (p)  |  [  — i -  +  — 


R  (Pfi,q)  R(P6,q)lT^  (p,q) 


<  a,62f  liLig^  -u(p)L  ds  ]2  + 

1  3D  RJ(P.,q) 


V2ff  H  dS]2 

j  3D  R(P  ,,  ,q)  R^_0t  (p  ,q) 


for  modified  constants  a^  and  &2 .  The  first  integral  on  the 
right  of  A-19,  with  Schwartz'  inequality,  becomes 


A-20 


_y(q)-u(p)  id  1  2 
3D  R3(p6,q)  q]  i 


ds 

3D  R  (Pfi,q)  J  3D  R  (P  ,q) 


Defining  Z, (p) ,  a  patch  on  3D  as. 


qe3D 

A-21  £d<P>  -  (q|  R(p,q)<d|) 


where  d  is  sufficiently  smaller  than  the  Lyapunov  radius  associated 
with  3D  so  that  A-9  and  A-10  may  be  used,  it  follows  that 

f  ^Srr  f  dS  r  dS 

A-22  -  q  +  q 

^D  R  (P^.q)  3D\Zd(p;  R^Tp^q)  '^d(p)  PTp^q) 


*  5* 


■d  r  2tt 

+  c  dp  £ - ,  3/2, 

o  Jo  (R  (P6,p)  +  o^) 


The  constants  A  and  c  are  of  no  particular  importance  for  our 
present  purpose  but  we  note  that  A  is  the  surface  area  of  3D. 

In  the  second  integral  in  A- 2 2  we  have  introduced  polar  coordinates 
(p,<j>)  on  the  disk  tangent  to  3D  at  p  and  used  the  inequalities 
(e.q.  Gunter  [4]) 


A-23  dS  <  2pdpd<{>  ,q£Z^(p) 


A-24  R(P .  ,q)  >  -/R2(P.,pl+R2(p,q)  >- 

P  -  2  6  ~  2 


R2(P6,P)+P2 


The  integral  in  A-22  is  easily  evaluated  and  we  have 


dS 

A-25  — rS 

3D  R  (P, 


,  8a  ,  r  i  i 

,q)  1  ?  L  1T^7ir  '  * 


R (Pgfp)  +d  ) 


1/2  < 


for  some  constant  c^.  Similarly  the  second  integral  in  (A-20) 
may  be  written 


a.„  f  lu^ruouv  <M  j3DNr  i»(<»j2^»<«i*  «« 

3D  R3(P5,q)  q  -  3DXZd(p)  d 


+  af 


u  (q)  -u  (p) 


JZd(p)  (92^+R2(p,q))3/2  q 


2  2  fd  '27T  0 !  u  ( p+q,) -U  (p) 

-  C1^U^L,  (3D)+  c2^w(p)  I  +  C3  dp  d*  77Z7T~T;  3/2 

2  o  o  5  +p  ) 


where  again  we  use  polar  coordinates  on  the  plane  targent  to  p  to 
parameterize  Z^fp)  and  employ  (A-23) ,  (A-24)  and  (3.21).  The 
vector  in  local  rectangular  coordinates  centered  at  p  is 

A-2?  =  (S,n,C(£,n))  ,  5=  pcos<i>  ,  n=  psin<|> 

where  ^  is  a  local  representation  of  the  surface,  hence  depends 
on  p,  but  also  satisfies,  since  d  is  smaller  than  the  Lyapunov 
radius 

,  ,  i+  a 

A-28  |  <ap  ,  a  =  Lyapunov  index. 


Hence  for  d  sufficiently  small 
>1-25  |  q1  i  <2p. 


The  second  integral  on  the  right  in  A-19  is  treated  similarly  yielding 


I  » 
.  vtj 
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d-32 


JiiMzHjgil  ds 

3D  R(P6,q)R2"C1(p/q)  qJ 


<bil  lul  l£20D)+b2lu(£>)  1 


fd  r 2tt  I  y  (p+q, )  -u  (p) 

+  b,  |  dp  d<J>  — < — - - 

3  J°  J“ 


Incorporating  A-30  and  A-31  into  A-19  we  have,  for  appropriately 
redefined  constants. 


A-32 


tli(q)-U<P)]  |n-[Y0(P6'q)"Y0(p'q)]dSq|2iCl61  lyl  iLoD) 
dD  q  ^ 


c2  6 lu (p) | 


+  c3  6 


d  217  p|  V  (p+q,) -y  (p)  j  2 
d<j>  ± 


dp 

o 


o 


c461+a  J  dp 


217  |u  (p+q,)-u(p)  |2 
d<|>  — - = - 

o 


Integrating  over  3D  and  employing  Fubini's  theorem,  which  applies 

1 


since  yeL, (3D)  and  both  — *— * — and  ,  /■»  -* — y- 

(e252+p2r/2  pi"Ve262+p2 


are 


integrable  in  the  p,4>  variables  for  6>0,  we  obtain 


A-33 


3D 


[y(q)-Vi(p)]|^r{Y0(P5,q)-Y0(p/q)  ]dSq|  dSp 


<  (c1+c2)6|  |u|  lIi2  (3D) 


+  c36 


rd  2+  p  rr  2 

!  p  d<f>  — 5-3 - j-ryv  |  |u(p+q1)-u(p)  I  dS 

>  0  Jp  (0  6  +P  )  • 
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+  c .  <5 
4 


i*»  (!  r?’  1 


[dp 


> 


d(f> 


x-ct 


o  'o  P  /e2<$2+p2 


>•  >  3D 


u (p+qx) -u (p) | 2dspJ 


The  integral 


3D 


2 

y(p+q^)-u(p) |  dSp  is  the  L2  modulus  of  continuity 
of  y  which  is  0(|q^|  )  (Stein-  [16  1 )  although  for  present  purposes 
it  is  sufficient  to  have  0(|q  |8)#  8>0.  Then  with  (A-29) 


A-34  |  | y (p+q. ) -u (p) |  dS  <  c  p 

3D  1  p 


3 


and  (A-33)  becomes 


A~2S  I  I  [  [u(q)-u(p)]|— [v  (p,,q)-y  (p,q)]dS  |2dS 
>  3D  J  3D  ~  °  Q  °  <1 


q1  p 


1  ^  «  Mull 


L2(3D) 


+  c2  6 


d 

p1+S  dp  +  c-51+ot 

o  (9252+p2) 3/2 


dp 

.l-a-0_ 


/P"62+p2 


%  . 
<  c  6 


2 

L, 


OD)+c2sS 


1+6 

7^f3/2 


ds 


1  o 


sl"°"S^2 


0^+s' 


When  3<1,  both  integrals  are  bounded  and  the  right  hand  side  is  0(<J8) 
If  S>1  even  more  rapid  decay  of  (A-34)  with  5  is  easily  shown.  In 
all  cases  the  validity  of  (A-2)  is  thus  ensured. 
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The  proof  of  A-3  involves  an  argument  similar  to  that  used  in 
establishing  A-2  which  we  merely  suggest.  With  (2.3)  (3.16),  (3.17) 

and  (A- 13)  we  obtain 


A-36 


■2-x(SyXPJ 

5nt 


3n~ 


(Sy)  (p)  r<4 


I  y  (q}  I  l3n“tY  (P<S  'q)  *Yo  (p6  'q) 
3D  p 
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The  first  integral  on  the  right  may  be  estimated  in  the  same  way 

3  3 

as  (A-17)  since  (A-16)  remains  valid  when  is  replaced  by 

q  p 

The  third  integral  is  precisely  the  same  as  (A-19)  and  is  treated 
identically.  The  second  integral  is  estimated  using  the  following 
inequality 
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which  enables  us  to  write 
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The  first  term  in  the  right  also  appears  in  A-19  and  is  treated 
as  it  was  there.  The  second  term  iso  (5°)  since  the  integral 

may  be  shown  to  be  0  (  .  This  completes  the  proof  of  Theorem  3.6. 

6 

We  remark  that  while  we  have  proven  the  theorem  for  Lyapunov 
surfaces  of  index  a,  the  proof  is  contingent  upon  the  existence 
of  families  of  parallel  surfaces,  which  we  have  only  asserted 
for  index  1.  This  is  sufficient  for  our  use  of  the  theorem  in 
Section  IV,  however  should  the  families  exist  under  the  weaker 
hypothesis  (see  Burago,  Maz'ja  and  Sapozhnikova  [1])  then  Theorem  3.6 
will  hold  in  that  case  as  well. 
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